Capillary wave dynamics on supported viscoelastic films: Single and double layers 
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We study the capillary wave dynamics of a single viscoelastic supported film and of a double 
layer of immiscible viscoelastic supported films. Using both simple scaling arguments and a con- 
tinuum hydrodynamic theory, we investigate the effects of viscoelasticity and interfacial slip on the 
relaxation dynamics of these capillary waves. Our results account for the recent observation of 
a wavelength-independent decay rate for capillary waves in a supported polystyrene/brominated 
polystyrene double layer [X. Hu et al, Phys. Rev. E 74, Of 0602 (R) (2006)]. 



I. INTRODUCTION 

Due to the presence of thermally excited capillary 
waves, the free surfaces of fluids Q|, complex fluids, 
and other generically soft materials are fluctuating struc- 
tures. Examining the surface dynamics of such soft ma- 
terials provides a window into their rheology. Conse- 
quently, capillary waves and interfacial dynamics have 
been probed by light scattering techniques 0, Q m a 
variety of systems, including membranes and monolay- 
ers 0, liquid metals [3, and polymer solutions 0, 0, 



brushes and gels [fj. These light scattering tech- 
niques probe surface dynamics at length scales in the 
micron range; however, newer experiments using x-ray 
photon correlation spectroscopy (XPCS) 0, ^| extend 
these measurements into the sub-micron range, allowing 
investigators to probe even smaller scale surface height 
fluctuations. 

The theory of the surface dynamics of polymeric ma- 
terials has focused on the continuum mechanics of vis- 
coelastic liquids, including two-fluid approaches to semi- 
infinite polymer solutions |12| and single component flu- 
ids supported on a rigid substrate 0, 0, These 
continuum approaches have been remarkably successful 
in accounting for the surface dynamics observed exper- 
imentally in many systems. In polymeric systems such 
continuum based methods must begin to fail in the limit 
of increasing molecular weight and decreasing film thick- 
ness, where the short dimension of the liquid becomes 
comparable to size of the constituent molecules. Indeed, 
there is experimental evidence ^|| for and subsequent 
theoretical conjecture 0,^3 about the shift (relative to 
its bulk value) in the glass transition temperature of high 
molecular weight polymer thin films. 

In recent experiments Lai, Sinha, and coworkers [l9l 
l2fj| used XPCS to study capillary waves on layers of 
polymeric films. They examined both single layers of 
polystyrene (PS) supported on a silicon substrate 0] 
and double layers of a PS layer overlying a brominated 
polystyrene (PBrS) layer on a silicon substrate [2(j|. In 
the double layer experiments, the XPCS technique al- 
lowed for the independent measurement of the dynamics 
of both the upper surface of the PS, which we refer to 
as the free surface, and the PS/PBrS interface, which we 



refer to as the buried interface. 

The most striking result of the Lai et al. double layer 
experiments is the appearance of a slow (~ 100 s) de- 
cay rate in the height-height correlations of both the 
upper and buried interfaces that is approximately inde- 
pendent of the in-plane wave number q. Such a phe- 
nomenon cannot be explained by previous theoretical 
work 0, 0, 0, 0| . In this article we extend the pre- 
vious theoretical work to two-layer viscoelastic systems 
and explore - via continuum visco-hydrodynamic calcu- 
lations and scaling arguments - the possible origins of the 
q-independent decay rate reported by Lai and collabora- 
tors. Acknowledging that the thickness of the fluid layers 
is on the order of a few radii of gyration of the constituent 
molecules, one may question whether these discrepancies 
point towards the failure of a continuum-based analysis. 
Despite this concern, we show that a viscoelastic contin- 
uum model can give rise to the observed q-independent 
decay rate. In order to quantitatively account for the ex- 
perimental data, however, we must use surprisingly long 
stress relaxation times in these continuum models, which 
in turn leads us to postulate that the polymer dynam- 
ics in the thin films is hindered by confinement effects. 
Based on these attempts to fit the scattering data using 
our model, we suggest that the experiments provide an 
interesting measure of confinement effects on the molec- 
ular dynamics in the melt. 

The remainder of the paper is organized as follows: 
in section [H] we examine simple scaling arguments that 
suggest that viscoelastic supported films can exhibit the 
phenomena reported by Lai and collaborators. In sec- 
tion lfffl we calculate the dynamics of the single supported 
fluid layer system using a continuum hydrodynamic de- 
scription of the fluid. We consider both a purely viscous 
Newtonian fluid and the simplest model for a viscoelas- 
tic fluid, the Maxwell fluid, which is characterized by a 
single relaxation time. We also examine the effects of 
slip at the fluid/substrate interface on these dynamics; 
in particular, we show that slip alone cannot account for 
the g-independent decay rate. We then turn to the study 
of the two layer system in section llVl Here, we consider 
both the system of two Newtonian fluids as well as the 
case of one Newtonian fluid and one Maxwell fluid. We 
show that our model of a Maxwell fluid buried beneath 
a Newtonian fluid can account for a number of the ex- 
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perimental features observed in the double layer systems 
studied by Lai et at. We again explore the effects of 
slip on the dynamics of the system, in this case at the 
liquid/liquid interface. Finally, in section Ivl we discuss 
our results more broadly in the context of the dynamics 
of multilayered systems, focusing on the long stress re- 
laxation times necessary to account for the results of the 
Lai group, as well as the implications such times have 
for the dynamics of polymers near an interface. We con- 
clude with suggestions for future experiments to further 
test our analysis. 
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FIG. 1: Schematic illustration of a single supported fluid 
layer of thickness d. 



II. SCALING ANALYSIS 



The experiments of Lai et al. 19, 20] present an inter- 
esting theoretical challenge for which at least a few sug- 
gestions have been offered 0, 0, Q| . We first use a few 
numerical estimates to narrow the focus of our problem. 
Consider a supported fluid film of thickness d on a solid 
substrate, as shown in Fig.^ The fluid has mass density 
p and a surface tension 7 at the free surface. The defor- 
mation of the free surface in capillary waves is subject to 
restoring forces due to gravity and surface tension. We 
assume that the bending energies of the interface are neg- 
ligible (we will return briefly to this point in section Ivjl . 
The importance of these two forces depends strongly on 
the wavelength of the dist urban ce. For wavelengths less 
than the capillary length y^f/pg, surface tension induced 
forces dominate over gravitational forces. The XPCS ex- 
periments can detect capillary waves of wavelengths less 
than the typical transverse coherence length of the beam, 

which is ~ io^m m m 113. The capillary wavelength, 

however, is on the order of 1 cm; therefore, we may ne- 
glect the effect of gravity. 

The importance of inertial effects in the fluid dy- 
namics may be estimated by considering the Reynolds 
number [2l| of the flows associated with the capillary 
waves. The Reynolds number for such flows is given by 
Re ~ Dihu/v, where D\ is the length scale over which 
the fluid velocity vanishes (we expect it to the be the 
lesser of the inverse wave number 1 jq and the film thick- 
ness d), v = rj/p is the kinematic viscosity and h, uj 
are the typical height and decay rate of the surface dis- 
turbance, respectively. Using the equipartition theorem 
to determine the average magnitude of thermally gener- 
ated capillary waves, (|ft-,j| 2 ) ~ ksT/^q 2 , we find that 
Re ~ 0(1O~ 7 ). This estimate demonstrates that inertial 
stresses are greatly dominated by viscous ones in the ma- 
terial. Thus, we ignore inertial stresses in the remainder 
of this article; that is, we assume that all of the fluids are 
completely overdamped. 

Given that we are considering low Reynolds number 
dynamics at scales well below the capillary length, we 
now turn to simple scaling arguments to determine the 
possible wave number dependence of the relaxation rates 
of overdamped capillary waves. In particular, we ask: 
what properties of the fluid lead to a ^-independent dis- 



persion relation, i.e. uj(q) ~ q°? The fluid deformations 
generated by capillary waves store energy in the fluid in- 
terface and, for viscoelastic fluids, in the bulk as well. 
When these deformations relax, this energy is dissipated 
through viscous stresses in the fluid. The scaling behav- 
ior of the dispersion relation can then be found by equat- 
ing the power generated by relaxing the elastic deforma- 
tions in the fluid with the power dissipated viscously by 
the fluid. 

For a Newtonian fluid, energy can be stored only in the 
interface at the free surface. Consider the lightly shaded 
portion of the fluid in Fig. ^ whose cross-sectional area 
in the x — y plane is A ~ l/q, I being a unit length in the y 
direction. The total power transferred to this volume of 
fluid from the interface is the product of the normal stress 
generated by the surface tension jq 2 h, the surface area 
A, and speed of the surface as the deformation relaxes, 
ujh: 



Psurf ~ {iq 2 h) Aujh. 



(1) 



The viscous force density for the incompressible fluid 
is 77V 2 v, where v is the fluid velocity. Then the power 
dissipated in the bulk of the fluid is given by 



d 6 xvr]V z v w Ar) I dzv(z)\7 2 v(z), (2) 




where v(z) is velocity averaged over one wavelength in 
the horizontal (x) direction; since the velocity is periodic 
in x, its variation in x does not affect the scaling behavior 
of the dispersion relation. 

We can relate the fluid velocity v to the interfa- 
cial height via volume conservation. During the re- 
laxation of a surface undulation, fluid volume is trans- 
ferred from elevated regions to the depressed ones, so 
Adth w J dydzv(z), where the integral is over a unit of 
area whose normal is parallel to x. Using this relation 
and equating the power input to the power dissipated, 
we find that the wave number dependence of the decay 
rate uj can be determined from 



11 



J dzv(z) 



77 J Q dzv(z)V 2 v(z) 



(3) 



3 



We now distinguish two different scaling regimes. In 
the thin layer limit (qd <C 1) the velocity decay into the 

fluid is set by the thickness of the layer, so dz ~ d and 
V 2 ~ l/d 2 . For thick layers (qd Cl), however, the ve- 
locity decays exponentially into the fluid, so J Q dz ~ 1/q 
and the V 2 ~ q 2 . Thus, we find two possible dispersion 
relations depending on the the value of qd, 




qd < 1 
qd > 1 



(4) 



This argument suggests that a Newtonian fluid will not 
exhibit the desired q-independent dispersion relation. 

If we consider instead a viscoelastic fluid, we note that 
elastic energy is also stored in the bulk due to the state 
of deformation in the material, as long as the decay rate 
of the capillary wave is fast compared to the stress re- 
laxation time in the viscoelastic material. This storage 
of elastic energy in the bulk changes our previous scaling 
arguments. For example, in an elastic solid character- 
ized by a bulk modulus /i, the power released in the bulk 
during the relaxation of the deformation takes the form 



P, 



bulk 



A\i \ v(z)V 2 u(z), 



(5) 



where u(z) is the average displacement field in the 
medium: u(z) ~ v(z)/u>. This expression is similar in 
form to the viscous power dissipation in Eq. (0). In- 
deed, if we equate the power dissipated with the power 
input from the bulk, we find uj ~ fi/rj for all values of 
q. Thus, this simple scaling argument suggests that the 
desired q-independent behavior is a manifestation of the 
viscoelastic response of the fluid. In order to verify the 
results of these heuristic arguments, and to determine 
over which range of wave numbers one might observe a 
q-independent dispersion relation, we now turn to the 
complete solutions of the Stokes equation for supported, 
overdamped Newtonian and viscoelastic fluid layers. 



III. SINGLE LAYER 

Consider the response of the single, supported, incom- 
pressible fluid layer shown in Fig. ^ to an external stress 
field a ext , which is normal to the free surface and char- 
acterized by frequency lo and wave vector q = gx, 



a e 



i(qx-iot) 



(6) 



The response of the fluid to the applied stress is described 
by the vertical deviation or height function of the free 
surface h(x,t), the velocity v(x,z,t), and the pressure 
P(x, z,t). Given the form of a ext in Eq. @), all of these 
dynamical quantities can be written in the form 



g(x,z,t)=g(q,z,u;)e i ^-" t \ 



(J) 



The general solution to the Stokes equation in the case 
that the dynamical quantities take the form of Eq. Q 



is determined in Appendix ^ In particular, all of the 
dynamical quantities listed above can be related to the 
normal component of the fluid velocity, which is given 
by Eq. i|A9() . In order to determine the four integration 
constants in Eq. (|A9|) . we need to specify the boundary 
conditions at the top and bottom boundaries of the fluid. 
At the fluid/substrate interface, the normal component 
of the fluid velocity must vanish, 



\z—d 



0. 



(8) 



We may account for both slip and stick boundary condi- 
tions on the tangential velocity component by introduc- 
ing a slip length A, so that at the fluid/substrate interface 



J *\z=d 



-\d z v x \ z=d . 



(9) 



Taking A = reduces the above boundary condition to 
the usual no-slip condition. 

At the free surface, the rate of change of the height of 
the interface must equal the fluid velocity at that point, 



-iujh(q, lj) 



(10) 



Furthermore, this interface cannot support shear stresses, 



f 



0, 



■xz\ z= o 

where the fluid stress tensor <jf takes the usual form 



(11) 



cr\, = rj(uj) (diVj + djVi) - SijP, 



(12) 



r)(u>) being the frequency dependent viscosity character- 
izing the viscoelastic response of the material. Finally, 
the free surface can support a stress discontinuity be- 
tween the externally applied stress and the hydrodynamic 
stresses in the bulk material, due to the presence of a fi- 
nite surface tension: 



~/q 2 h(q,uj) 



J f I 



(13) 



Using the boundary conditions Eqs. (|8|)- (|llfl to solve 
for the fluid velocity and pressure fields, we obtain the 
normal fluid stress at the free surface from Eq. I|12l) . 



where 

B(q) = Aq 



a L\ z=0 = iur)(u)B(q)h(q,u), (14) 
cosh 2 (qd) + q 2 d 2 + qX (2qd + sinh(2gd)) 



sinh(2qd) - 2qd + iq\ sinh (qd) 



Then Eq. (jT3*j) becomes 



(T = h(q,u) [jq 2 - iujf](uj)B(q)] . 



(15) 



(16) 



The normal mode frequencies are those which satisfy 
Eq. I |lti|l in the absence of an external stress and with a 
non-zero fluid height: 



7<? 2 - iuj n r)(u n )B(q) = 0. 



(17) 
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FIG. 2: Dimensionless decay rate rjdui/'y for a Newtonian fluid 
as a function of qd. The solid (dotted) lines are for a system 
with a small (large) liquid/substrate slip length, A = O.Old 
(A = Wd). 



For this overdamped system the frequencies u> n , given by 
the solutions to Eq. I|17(l . all lie on the negative imaginary 
axis. We refer to the norm of these complex numbers as 
the decay rates of the system, Co n — iu> n . 

The experimentally measurable quantity for this sys- 
tem is the intensity autocorrelation function, which is di- 
rectly proportional to the height-height correlation func- 
tion S(q,t) oc (h(q,t)h*(q,0)) (see Appendix FBI for fur- 
ther details). Using Eq. I|B3|) . we find the height-height 
correlation function to be 



S(q,t) = - V lim 



(w - uj„ ) er 



to [jq 2 — iujr](uj)B(q)] 



(18) 



The decay rates in the presence of both small and large 
slip lengths are shown in Fig. [21 When the slip length 
is small, A -C d, the decay rate is independent of it to 
leading order in q, 




qd < 1 
qd > 1 



(20) 



Other than numerical prefactors, this result is identical 
to the one obtained using the simple scaling arguments 
given above - see Eq. J2J. It is also consistent with the 
experimental results obtained by Sinha et. al for a single 
PS layer, for which a q 4 scaling behavior was observed 
over wave numbers qd < 1 [l9l | 

When a significant amount of slip occurs between the 
fluid and the solid interface (i.e. A 3> d), an intermediate 
scaling regime in the decay rate appears, 



ui(q) 



ri£jL2 qd<l, q 2 d\<^\ 
^ gd«l, 9 2 dA»l 



(21) 



UL 

2n 



qd > 1 



Thus, we can see that a Newtonian fluid does not exhibit 
a g-independent decay rate, even if there is a significant 
amount of slip between the fluid and substrate. 

Because there is only one decay rate, the height-height 
correlation function Eq. I|18|) exhibits a simple exponen- 
tial decay, 



S(q,t) 



1 

iq 2 



(22) 



Thus, the surface height dynamics - as parameterized 
by the decay rates uj n - can be extracted from the ex- 
perimental measurement of S(q,t). We now turn to the 
calculation of these decay rates for single supported New- 
tonian and Maxwell fluids. In light of the Lai et al. data 
we pay particular attention to the dependence of these 
decay rates on wave number q. 



A. Newtonian Fluid 

For Newtonian fluid the viscosity is real, positive, and 
independent of frequency, r/(uj) = 77 > 0. This is gen- 
erally a good approximation for small molecule liquids 
and for viscoelastic materials on time scales much longer 
than their typical stress relaxation times. As shown in 
Appendix El we expect only one decay rate, and indeed 
there is only one root of Eq. (|17fl in this case, 



u>(q) 



iq 

T)B(q)- 



(19) 



We note that, for the case where there is no slip between 
the fluid and substrate, the same result can be obtained 
from the previous theoretical work of Jackie |13l Il9| . 



B. Maxwell Fluid 

A Maxwell fluid, which is the simplest model for a vis- 
coelastic material, has a complex, frequency dependent 
viscosity of the form 



Et 



1 — iujt ' 



(23) 



where E is the transient modulus of the polymer network, 
r is the stress relaxation time of the medium, and 77 ac- 
counts for the high frequency viscous response p^ . In 
this case Eq. <|17[) can be written as a quadratic equation 
in oj; that is, we obtain two decay rates, in agreement 
with the arguments given in Appendix El These rates 
are given by 



w±(g) 



A/~± \JN 2 - 4f3r]T 
2^7 ' 



(24) 



where N = r\ + Et + (3t and (3(q) = r yq 2 /B(q). Since 
Af 2 - 4f3r]T = (77 + Et - (3t) 2 + 4/3 Et 2 > both the fast 
decay rate a) + and slow decay rate cD_ are positive real 
numbers. 

In Fig. |31 we plot the fast and slow decay rates for both 
small and large amounts of slip at the substrate, using 
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Parameter 


symbol 


Value 


PS layer thickness 


d 2 


100 nm 


PBrS layer thickness 


d [SL], di [DL] 


200 nm 


PS viscosity 




10 4 kg/m-s 


PBrS viscosity 


V [SL], m PL] 


10 6 kg/m-s 


PBrS surface tension [SL] 


7 


10" 2 N/m 


Interface tension PS/PBrS 


7i 


10~ 3 N/m 


PS surface tension 


72 


10" 2 N/m 


PBrS Plateau modulus 


E 


10 3 Pa 


PBrS stress relaxation time 


T 


100 s 



TABLE I: Geometric and rheological parameters correspond- 
ing to the single layer [SL] and double layer [DL] Lai exper- 
iments |20|. Both the plateau modulus and the stress relax- 
ation time are chosen so as to reproduce the major features 
observed experimentally in the double layer system (see sec- 
tion E). 



parameter values characteristic of a layer of PBrS, as 
listed in Tabled When the slip length is small compared 
to the film thickness, A <C d, we find 



{ V+Et 
2)7 



<D_(g) = 



= I 3(i,+Et) 



qd < 1 
qd^> 1 ' 

qd < 1 
qd > 1 



(25) 



(26) 



When the slip length is large compared to film thickness, 
A 3> d, we find 



w_0?) 




7 g 4 rf 2 A 
■q+Er 



4(?H 
1 



Et 



< 1, 9 2 rfA < 1 
qd < 1, q 2 dA > 1 
qd > 1 



(27) 



(28) 



As we can see from both the analytic expressions and 
the numerical results shown in Fig. the fast (slow) de- 
cay rate is independent of q at low (high) wave numbers. 
Such behavior is in agreement with the experimental re- 
sults of Lai et ai, and it is also consistent with the scaling 
arguments given in section [n] for the elastic response of a 
viscoelastic fluid. On the other hand, the scaling of the 
fast and slow decay rates in their respective (/-dependent 
regimes - at high qd in the former case, and at low qd in 
the latter - is identical to that of a Newtonian fluid in 
these regimes. Thus, the viscous response of the fluid de- 
termines the scaling of Cu± where each decay rate exhibits 
a strong dependence on the wavelength, whereas the elas- 
tic response of the fluid gives rise to the g-independent 
behavior of lu± . 




FIG. 3: Dimensionless decay rates r]du±/y for a Maxwell fluid 
as a function of qd, using the parameter values listed in Table|l] 
for PBrS. The solid (dotted) lines are for a system with a small 
(large) liquid/substrate slip length, A = O.Old (A = 10d). The 
dashed lines indicate the approximate window of decay rates 
that can be measured in the experiments, 10 _3 s _1 < Co < 
10s" 1 . Inset: Ratio of the amplitude of the term in S(q,t) 
that decays with rate uj+ to that of the term that decays with 
rate Hi-. 



From this analysis we predict that the height-height 
correlation function S(q, t) will exhibit a double expo- 
nential decay for a Maxwell fluid. In particular, Eq. <B3I) 
becomes 



(W±T - 1) 



ryrB(q)Lb ± (Q± - iu T ) 



(29) 



These amplitudes are plotted in the inset of Fig. [3J 

In order to predict the form of experimental measure- 
ments of S(q, t), one must consider the values of both the 
decay rates and the (/-dependent amplitudes appearing in 
Eq. (|29|l . In principle, the decay of the height-height cor- 
relation function S(q, t) should always have the double 
exponential form given in Eq. It is clear fr om the 

inset of Fig. |3 however, that the amplitude of one of these 
terms can dominate the other at certain wave numbers. 
In particular, the amplitude of the slow mode dominates 
for qd <^1, whereas the amplitude of the fast mode dom- 
inates for qd ^> 1. In these regimes, it will be hard to 
measure both decay rates. The full double exponential 
decay should be observable only in the crossover region 
between these two regimes, which occurs at qd ~ 1 for 
parameter values consistent with current experiments. 

We can see from the inset of Fig.|3|that for both qd <C 1 
and gd > 1, the amplitude of the q-dependent rate is 
much larger than that of the q-independent rate. Taken 
in isolation, this fact suggests that it would be difficult to 
observe the ^-independent rate in these regions. In pre- 
dicting the experimentally observed behavior of S(q,t), 
however, it is important to take the value of the decay 
rates themselves into account. In particular, there is a 
finite range of decay rates that can be measured: the 
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h 2 (x,t) 



z=0 



z=d~ 



z=d 




FIG. 4: Schematic illustration of two fluid layers, of thick- 
nesses di and 4, on a substrate. The coordinates are chosen 
so that the fluid/solid interface is the plane z — d and the 
unperturbed free surface is the plane z = 0. 



decay times can be too fast or too slow to be observed 
within the time scales of the experiment. We have in- 
dicated a reasonable experimental window - from 1/10 
sec to 1000 sec for the decay times - by the dashed lines 
in Fig [21 Then we see, for example, that the slow decay 
rate becomes too slow to be measured experimentally at 
long wavelengths. As a result, the faster q-independent 
decay rate could be observed in this regime, despite the 
fact that its amplitude is much smaller than the slow de- 
cay rate. We do not comment on the sensitivity of the 
experiments to small amplitude capillary dynamics. If 
the larger amplitude mode is too slow and the smaller 
amplitude mode generates too small a surface height un- 
dulation, no interfacial dynamics may be detected. 

This analysis suggests that for qd < 1 it should in 
principle be possible to observe a q-independent decay of 
S(q, t) over at least one decade of qd, depending on the 
sensitivity of the measurement. For higher wave num- 
bers, qd > 1, however, we expect this g-independent be- 
havior to be obscured by the slower q-dependent mode 
with a larger amplitude, since it will then be fast enough 
to be measured by the current experiments. 



IV. DOUBLE LAYER 

We now turn to the double layer case, illustrated in 
Fig. 01 The buried and upper fluids have viscosities 
771(0;), 772(0;) and thicknesses d\,d 2 , respectively. The to- 
tal thickness is d = d\ + d 2 . In general, we can drive 
this system by spatially oscillatory normal stresses act- 
ing on both the buried fluid/fluid interface and the free 
surface. As in the single layer case, we consider stresses 
of the form of Eq. @ . Combining the stresses on both 
interfaces into a single vector, we write 



W) 



J(qx-ut) 



(30) 



For simplicity, we consider the possibility of slip at only 
the fluid/fluid interface, as the effects of slip at the 



fluid/substrate interface have already been explored in 
the single layer case. 

At the fluid/substrate interface all components of the 
fluid velocity must vanish: 



«i,*U rf = 0. 



(31) 



At the fluid/fluid interface the ^-components of the ve- 
locities of the two fluids must match, but the tangential 
components of the velocities do not if there is a finite slip 
length A at this interface: 

vi.z\ z=d2 = u 2, z L =d2 , (32) 
(v 2 ,x - vi,x)\ z=d2 = -A (d z v 2 , x ~ d z v hx )\ z=d2 .(33) 

Note that A > when the buried polymer layer is more 
viscous than the upper layer (Re [771(0;)] > Re [772(0;)] for 
all real ui), and vice- versa. 

From the definition of the interface heights, 



V2,z\ z=0 = -iuh 2 (q,uj). 



(34) 
(35) 



In addition, the free surface cannot support shear 
stresses, and the shear stress must be continuous across 
the fluid/fluid interface: 



7 2,xz 



z=0 



0. 



= 4,xz 



(36) 
(37) 



Finally, the normal stress discontinuities at both the 
fluid/fluid interface and the free surface are determined 
by their respective surface tensions (71 and 72, respec- 
tively), 



0-1 
02 



7ig 2 /ii(g,cj) 



z=d 2 



r )2q 2 h 2 {q,u)) - <r[ 



2,zz 



z=0 



(38) 
(39) 



Using Eqs. (|31|l - (|37|l to eliminate the fluid velocities in 
favor of the interfacial height functions, the normal stress 
equations Eqs. (|38|l and (|39|l can be written in matrix 
form: 



<x = S • h, 



(40) 



where cr 



£ = 



h 2 (q,uj) 



■ 2 _ ^sn(g^) 
' iq A(g,w) 

A(q,u) ^2<7 



, and 



A(q,w) 
2 IU)S 2 2 (<?.") 



(41) 



with 
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s n (q,w) =Aqrjl [l + (l + 2q 2 d\) sech (2qd 1 )] [tanh(2gd 2 ) - 2qd 2 sech (2qd 2 )} + Mvl [l - (l + 2q 2 dj) sech(2gdi)] 
* [tanh (2qd 2 ) + 2qd 2 sech (2qd 2 )] + 8qrjirj 2 [2qdi (2q 2 d\d 2 — l) sech {2qd\) sech (2qd 2 ) + tanh (2gdi)] 

+ 8g 2 A (77a - r? 2 ) {r?i [l + (l + 2q 2 d\) sech (2gdi)] [1 - sech {2qd 2 )\ (42) 
+ rj 2 [tanh (2qdi) — 2qd\ sech (2gdi)] [tanh (2qd 2 ) + 2qd 2 sech (2gd 2 )] |, 



Si2(g,w) = - 



2cosh(<7<i 2 ) ^ sech ((7(^2) 



{772 [l - (l + 2q 2 d 2 ) sech (2qdi)] [qd 2 + tanh (qd 2 ) 



2q i d{d 2 sech (2gdi) + [tanh (2qdi) - 2qdi sech (2qdi)] [qd 2 tanh (gd 2 ) + 1] 



2q\ (rji — rj 2 ) [tanh(2gdi) — 2qd\ sech(2gdi)] [tanh (qd 2 ) + qd 2 ] j, 



(43) 



s 22 (q, u) =4qr)2 t 1 _ i 1 + 2 <7 2 d 2 ) sech (2qdi)] [tanh (2gd 2 ) + 2qd 2 sec (2gd 2 )] 

+ 4g?7i?7 2 [tanh (2gdi) - 2gdi sech {2qd\)\ [l + (l + 2q 2 d\) sech (2qd 2 )] 

+ 8q 2 Xr] 2 (771 — J72) [tanh (2gdi) — 2gdi sech (2gdi)] [tanh (2qd 2 ) + 2qd 2 sech (2qd 2 )] 



(44) 



A(g, w) =2?72 [l - (l + 2q 2 d\) sech (2gdi)] [1 - sech (2gd 2 )] + 27/i [tanh (2gd x ) - 2qdi sech (2gdi)] (45) 
* [tanh (2qd 2 ) — 2qd 2 sech (2gd 2 )] + 4q\ (7/1 — i] 2 ) [tanh (2gdi) — 2qdi sech {2qd\)\ [1 — sech (2qd 2 )\ , 



where, for clarity, we have suppressed the (possible) fre- 
quency dependence of the viscosities 771 and rj 2 . 

The two normal modes of the double layer system are 
easily identified if we diagonalize the normal stress matrix 
equation Eq. (|41|l . bringing the dynamical relations into 
the form 



(46) 



0+ 




A+ 




h+(q,ui) 






A_ 




h-(q,ui) 



h±=A±h 1 +A 2 L h 2 , (47) 



where 

a± = Afai + A^o- 2 , 

and A± and A^ are, respectively, the eigenvalues and the 
orthonormal eigenvectors of the matrix S. The eigenval- 
ues may be written as 



Tr S ± y (Tr £) 2 -4det£ 



(48) 



and the corresponding eigenvectors are 



A± = 



1 



^12 



(A±-E u ) 



£i 2 
A± - Sn 



(49) 



We can see from Eq. (|46|) that h+ and h- are the am- 
plitudes of the two independent normal modes of this 
double layer system. As in the single layer case, the char- 
acteristic decay rates of these two modes can be found by 
solving the system of equations in the absence of external 
forces, cr = 0. From Eq. I|46|) we note that the character- 
istic decay rates Co^ of the h± mode are the roots of the 
eigenvalue A±, where n indexes the roots. 

Examining Eq. I|48|l , it is clear that these roots are also 
roots of the determinant of S, which may be written as 



det £(w) 



$(g,<") 
A(q,u) 



M 2 



A(g,w 



[7i s 22(<7, w) + J 2 su(q,uj)} ioj + 7i7 2 q 4 



(50) 
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where 

<£>(<?, w) =16g 2 ?7i?72 [4g 2 did 2 (l - q 2 d 1 d 2 ) sech(2gdi) sech(2<7d 2 ) + tanh (2qd 1 ) tanh (2qd 2 )] 
+ &q 2 ril [l - (1 + 2q 2 dj) sech(2qd 1 )] [l - (l + 2g 2 df.) sech(2 9 d 2 )] 
+ 8q 2 Vii£ [!+(! + V^i) sech(2gdi)] [l + (l + 2g 2 dl) sech(2gd 2 )] 

+ 16g 3 Ar/ 2 - 772) |r?i [tanh (2 g d 2 ) + 2 9 d 2 scch (2gd 2 )] [l + (l + 2q 2 d\) scch (2gd 2 )] (51) 
+ r/ 2 [tanh (2gdi) - 2qd x sech (2gdi)] [l - (l + 2<7 2 d 2 ) sech (2gd 2 )] } 



Furthermore, we can see from Eq. (|48|l that the trace of 
S is positive (negative) for the roots H>~ (o) + ). 

Given the normal modes of the system, we can again 
use the results of Appendix [B] to calculate the height- 
height correlation functions S%(q,t) oc (hk(q, t)ti£.(q, 0)) 
for k — 1,2, where k = 1 labels the buried interface and 
k = 2 the free surface. Eq. (|B3|> gives the normal mode 
correlation functions S± (q, t) , which can then be used to 
compute Sk{q,t) using Eq. (|4"7|l . 



S fc ( g ,t)=-£ 



la;,. 



d(ioj) 



(det S) 



where fc = 2, 1 for fc = 1, 2. 

Clearly for the two layer system our solutions depend 
on a larger set of geometric and material parameters. We 
do not show all possible parameter regimes, but consider 
in both the following figures and asymptotic results a 
parameter regime consistent with the experiments of Lai 
et al. 20] . In particular, we take the layer depths to 
be of the same order, di ~ d 2 , and the surface tension 
of the fluid/fluid interface to be less than that of the 
free surface, 71 < 72. We also assume that the buried 
polymer layer is much more viscous than the upper layer, 
Re [771(0;)] 3> Re[?7 2 (aj)] for all real to. However, many 
of our results below (e.g. the scaling behavior of the 
normal mode decay rates) apply more generally to two- 
layer viscoelastic systems. 



Two Newtonian Fluids 



In the case that both layers can be described as New- 
tonian fluids with viscosities 771(0;) = 771 and T7 2 (o;) = ?7 2 , 
it is clear from Eq. (|5Ull that det S is a simple quadratic 
function of to. Thus, there are two normal mode de- 
cay rates, as expected from the arguments given in Ap- 
pendix [C] Specifically, the roots of Eq. (|50f) are given 



— = — CO 




FIG. 5: Dimensionless decay rates r]\duj ± /71 and amplitude 



(52) ratios A + /A + (inset) as a function of qd for a double layer 
system with a buried Newtonian PBrS layer and an upper 
Newtonian PS layer, using the parameter values given in Ta- 
ble U The solid (dotted) lines are for a system with a small 
(large) liquid/liquid slip length, A = O.Old (A = lOd). 



by 



£5 ± («) 



q_ 

2$ 

± 



71S22 + 72S11 



(53) 



V (7l s 22 



72S11) — 47i7 2 $A 



where tD* = 



; ± . The values of the fast mode ui + and 



the slow mode u> for parameters consistent with a PBrS 
buried fluid and a PS upper fluid (see Table[U are plotted 



in Fig. The u> 



r,+ 



modes correspond to nearly in 



phase and out of phase undulations of the two surfaces, 
respectively. 

Using Eqs. (|42|I - H45|) . one can show that both normal 
mode decay rates exhibit the same asymptotic scaling 
behavior observed for the single viscous layer discussed 
above. When the slip length is small, A <C d, 



-(g) 



"(9) 



72 q d 2 

3»72 

72k! 

2))2 




qd < 1 
qd > 1 

qd < 1 
qd > 1 



(54) 



(55) 
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A large slip length at the fluid fluid interface (i.e. A ^> d) 
produces an intermediate scaling regime in uj + , 

^ + (l) = \ 22 ^ qd<l,q 2 dX>l. (56) 

However, the value of Q~ is unaffected by the presence 
of a large slip length when r]i ;» r/ 2 . Thus, we can see 
that, like the single layer double layer of Newto- 

nian fluids does not exhibit a g-independent decay rate, 
even in the presence of liquid-liquid slip. 

Examining the normal mode amplitude ratio as a func- 
tion of wave number, shown in the inset of Fig. El we see 
that at the upper surface the fast mode amplitude al- 
ways dominates, while at the buried interface the slow 
mode always dominates. Thus, the decay of the height- 
height correlation functions, which are given by Eq. (|52H . 
are essentially of the form of a single exponential decay, 
Si(q,t) ~ e _£rt and S 2 (<?,*) ~ e~" +t . The large size of 
the amplitude ratios at the two surfaces is due primarily 
to the large viscosity difference between the two layers. 
If we were to consider a two layer system having similar 
viscosities and similar interfacial tensions at the two sur- 
faces, we would find amplitude ratios at the two surfaces 
of order unity. Of course, the rates themselves would be 
approximately equal as well in this limit. 



B. One Maxwell Fluid, One Newtonian Fluid 

When one of the fluids in the double layer is viscoelas- 
tic, we expect to find additional decay rates, as we found 
for the single viscoelastic layer. Here, we focus on the 
case of a buried viscoelastic fluid. It is known from the 
single layer experiments that PS fluid layers - the upper 
fluid in the double layer experiments - are well described 
by a purely viscous Newtonian fluid model [l^. On the 
other hand, we expect that the PBrS layer to have longer 
stress relaxation times. We model this viscoelastic layer 
using the Maxwell model, so that 771(0;) is of the form of 
Eq. This approach introduces two unknown rheo- 

logical parameters for the material (namely, E and t), 
which we discuss more fully below. For simplicity, we 
analyze this system with stick boundary conditions at all 
interfaces, i.e. A = 0. 

From examining the zeros of Eq. (|50|l we find four sep- 
arate decay rates for this dynamical system, in agreement 
with the arguments given in Appendix El These decay 
rates for a double layer system with a buried PBrS layer 
and an upper PS layer are shown in Figs. EK and EP: 
using the parameters given in Table [i] (these two plots 
are identical, except that the line styles in each are cho- 
sen to correspond to Figs. HJ3 and |Sf>, respectively; we 
explain this in more detail below). The behavior of the 
decay rates is qualitatively similar to the single layer vis- 
coelastic case. In particular, each rate has two major 



scaling regions: a g-dependent region, where the decay 
rate scales like the corresponding one in the Newtonian 
fluids case (i.e. ~ g 4 for small qd or ~ q for large qd), 
and a g-independent region. In the single layer case, how- 
ever, the crossover between these two regions is abrupt 
and occurs at qd ~ O(l). This is not the case for the 
double layer system; indeed, one decay rate is essentially 
constant over the entire range of wave numbers shown in 
Fig. El with its g-dependent regime not appearing until 
higher wave numbers (data not shown). 

In principle, all four decay rates shown in Fig. Eicon- 
tribute to the height-height correlation functions at both 
the upper and buried interface. As in the single layer 
case, though, the amplitudes corresponding to each de- 
cay rate are important in determining which rates can 
be observed experimentally in these correlation functions. 
Figs.EjB and Ef> show the amplitudes for each decay rate 
at the buried and upper interfaces, respectively. The de- 
cay rates shown in Figs. EJV and EP are plotted using 
the same line style (dotted, black, blue, and red) as their 
corresponding amplitude at the buried and upper inter- 
faces, respectively. In these two figures we represent the 
experimentally accessible range of timescales by the re- 
gion between the two grey dashed lines. From examining 
Figs. EJ3 and EP it is clear that the ratios of the vari- 
ous amplitudes vary by many orders of magnitude over 
the experimentally accessible wave number range. Also, 
the dominant amplitude at both interfaces is always q- 
dependent, as was seen in the single layer case. 

We consider first the dynamics of buried interface 
shown in Figs. EJA andEP- If we focus on the wave num- 
ber regime accessible in the experiments of Lai et al, 
0.8 < qd < 2, we see that the slow, g-dependent rate, 
which has the largest amplitude at the buried interface, is 
too slow to be detected experimentally. The other decay 
rates all have amplitudes of approximately the same mag- 
nitude in this region. Thus, it appears that the buried 
layer dynamics should, in fact, be best described by a 
height autocorrelation function having multiple exponen- 
tial decays. The data of Lai et. al, however, was con- 
sistent with single exponential decays at this interface. 
This could be due to the fact that decay rates do not dif- 
fer substantially in this region (see the blue, black, and 
dotted curves in Fig. EK)i thus rendering it difficult to 
observe the multi-exponential decay profiles. We do pre- 
dict, however, that the g-independent decay rate should 
be observable at the buried layer, in agreement with the 
experimental data. 

We now turn to the dynamics of the free surface 
shown in Figs. EP andEp. We see that for wave num- 
bers qd ~ 0.7, it should be possible to observe a dou- 
ble exponential decay, since the amplitudes of the dom- 
inant, g-dependent decay rate (red line) and subdom- 
inant g-independent rate (blue line) arc approximately 
equal, as shown in Fig.|f)p. Furthermore, we can see from 
Fig. EP that the g-dependent rate is slower than the g- 
independent rate in the experimental range of wave num- 
bers, 0.8 < qd < 2. The experiments do indeed observe a 
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FIG. 6: {color online) Dimensionless decay rates rjid^/ji (A, C) and amplitudes ^iA/d 2 at the buried (B) and upper (D) 
interfaces, as a function of qd for a double layer system with a buried Maxwell PBrS layer and an upper Newtonian PS layer in 
the absence of slip, using the parameter values given in TableU] (^4, C) The rates shown in these two figures are identical, except 
that the line styles in Fig. A correspond to those in Fig. B, while the line styles in Fig. C correspond to those in Fig. D. The 
dashed lines indicate the approximate window of decay rates that can be measured in the experiments, 10 _3 s _1 < Co < 10s _1 . 
(B, D) The amplitudes in these figures are indicated by dotted, black, blue, and red lines, in order of increasing amplitude. 



double exponential decay with a faster g-dependent rate 
and a slower g-independent rate; however, the amplitude 
of the g-independent decay is larger than that of the q- 
dependent decay in the experiments, whereas the oppo- 
site is true for our theoretical predictions. 

We suspect that pursuing more detailed numerical fits 
of the theory to the current experiments is unproductive 
due primarily to the inadequacies of the Maxwell model. 
In particular, a more realistic description of the stress 
relaxation in the PBrS layer will allow for a spectrum of 
relaxation times, leading to a band of decay rates at a 
given wave number. As is clear from Figs.lHP and D, the 
amplitudes of the decay rates are highly nonlinear, and 
thus are sensitive to the details of the viscoelastic model 
used for the fluid. As the theory currently stands, it ap- 
pears to be fortuitous that one may observe the double 
exponential decay at the free surface, since amplitude 
ratio between the two dominant modes (red and blue 
lines) in this figure have an extremely narrow maximum 
at qd k 0.7. We expect, however, that a broadening 
of the relaxation time spectrum in the material will in- 
crease the width of the features in Fig. |SjD, making the 
observation of the multiple exponential decay in the data 
much more plausible. It is important to note that the use 
of a more realistic viscoelastic model for the buried fluid 



should not affect the existence of the g-independent decay 
rates themselves. Indeed, the scaling arguments given in 
section [H] show that the presence of the g-independent 
decay rate is a robust feature of the elastic response of 
the fluid. Thus, a more realistic model for the buried 
fluid should be able to account for the discrepancies be- 
tween our theory and the experimental results of Lai et 
al. 

Finally, we note that when the upper fluid is a Maxwell 
fluid and the lower fluid is a Newtonian fluid, an addi- 
tional decay rate appears, in agreement with the argu- 
ments given in Appendix El Although four of the rates 
behave in a manner similar to that of the four decay rates 
for a buried Maxwell fluid - with each rate displaying one 
region of g-independent behavior and another region of 
g-dependent, Newtonian-like behavior - the fifth rate is 
essentially constant over the wave numbers of interest. 
Furthermore, we find that at certain wave numbers the 
amplitude of one of the g-independent decay rates can ac- 
tually be as large as the largest g-dependent amplitudes 
in the height-height correlation functions, at both the 
buried and upper interfaces (data not shown). Thus, our 
analysis suggests that detecting the g-independent rates 
may be even easier when the upper fluid is viscoelastic. It 
is important to repeat, however, that we do not suspect 
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that it is a viscoelastic response of the upper PS layer 
that is leading to the ^-independent decay rates seen in 
the experiments, since the purely viscous model fits the 
single layer PS data well Qjj. 



V. CONCLUSION 

In this article we have examined in detail the over- 
damped dynamics of purely viscous and viscoelastic sup- 
ported films, taking the effects of both liquid/substrate 
and liquid/liquid slip into account. These calculations 
show that a simple approach to understanding the ap- 
pearance of wave number independent decay rates in the 
height autocorrelation function S(q, t) in both single and 
double layer films is found by allowing for a viscoelastic 
response of the material. Simple scaling arguments show 
that a viscoelastic material will exhibit such dynamics 
over some finite wave number range. The more detailed 
continuum mechanics calculations based on the Maxwell 
fluid model presented above show that this wave number 
independent scaling regime is experimentally accessible, 
at least for a viscoelastic material with the appropriate 
values of the plateau modulus E and the stress relaxation 
time t. 

More generally, the appearance of a ^-independent re- 
laxation rate in S(q, t) appears to be a signal of a vis- 
coelastic response of the supported polymer film on the 
time scales probed by the measurement. The observed 
decay time is strongly controlled by the longest stress 
relaxation time in the material, and thus may serve as 
an important measure of the rheological properties of 
such thin supported films. In the numerical results that 
we presented our choice of the stress relaxation time r 
is constrained by the value of the g-independent decay 
time measured in the experiments, which is of the or- 
der of 100s. In order to fit the data, we had to take r 
to be of this same order, which is almost two orders of 
magnitude larger longer than might be expected for PS 
at the reported molecular weight [23j, We suggest that 
discrepancy may be attributed to two factors. First, the 
bromination of the PS changes the microscopic dynamics 
of dynamics of the polymer. Second, the narrow thick- 
ness of the PBrS layer (the layer is about three radii of 
gyration of the chain in ^-conditions) adds further con- 
straints to the chain dynamics leading to a longer stress 
relaxation time. 

Given a determination of r, we found that our choice 
of the plateau modulus was not as precisely constrained 
by the data. If we choose Er/r/i 3> 1, the value of the 
g-independent decay time is decreased below r. Since the 
value of r needed to fit the data is already anomalously 
large, we wish to choose the smallest possible value of 
t: this restricts us to the region Er/rji C 1. As long 
as we choose a value of E within this region, though, 
the value of the g-independent decay rate is unaffected 
by it. Rather, the principal effect of varying the plateau 
modulus is to shift the various amplitude ratios of the 



multiple decay rates, in particular the wave numbers at 
which these ratios approach unity. 

This points out two more general results. The first is 
that the appearance of multiple exponential decay rates 
of the height autocorrelation function S(q, t) for a sin- 
gle layer system is a generic feature of our viscoelastic 
model. In a two-layer system, even purely viscous fluids 
admit a double exponential decay. If either of the layers 
is viscoelastic, then this double decay becomes a more 
complicated multiple decay as discussed above. Based 
on these calculations we suspect that the observation of 
multiple exponential decays of S(q, t) in a single layer 
system may be taken as evidence of the viscoelasticity of 
the supported film. Secondly, our calculations show that 
observable multiple exponential decays depend on am- 
plitude ratios that are in turn controlled by the plateau 
modulus of the material. These ratios can vary by many 
orders of magnitude with the plateau modulus and wave 
number. The observation of multiple exponential decays 
requires the ratio of the two largest amplitudes to be 
near unity; from our Maxwell-model based calculations 
such occurrences occupy a small part of the phase space 
spanned by wave number and plateau modulus. While 
we expect that a broader spectrum of stress relaxation 
times in the material will enlarge the region of phase 
space over which these multiple exponential decays can 
be seen, such multiple exponential decays may not be a 
generically observable feature of viscoelastic supported 
films. When such multiple exponential decays are ob- 
served, however, they should provide a sensitive window 
onto the plateau modulus of the material and thus mea- 
sure the entanglement length in the layer. 

Finally, we point out that these calculations do not in 
general exclude other potential mechanisms as the under- 
lying cause of the dynamics as observed by XPCS. They 
do show, however, that these data are not consistent with 
the dynamics of two immiscible Newtonian fluids either 
with stick or slip boundary conditions at their bound- 
aries. Moreover, one can show that postulating more 
complex surface energy functionals - including, for ex- 
ample, an interfacial bending modulus - will only lead 
to relaxation rates of S(q, t) that are even more strongly 
dependent on the wave number, which would be incon- 
sistent with the data. Therefore, such considerations can 
also be excluded. It appears that most simple way to ac- 
count for the data is to postulate a viscoelastic response 
of the supported film with a stress relaxation time on the 
order of the observed decay rate of the height autocorre- 
lation function. 

A more detailed analysis of the interfacial dynamics 
of complex fluids via XPCS holds the promise of prob- 
ing molecular motion in confined geometries that may be 
interpreted with the aid of the continuum modeling pre- 
sented in this article. If our calculations are to serve in 
this manner, however, they must first be tested further 
by experiment on rheologically well-characterized mate- 
rials of thickness large enough to discount the effects of 
molecular confinement. 
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APPENDIX A: SOLUTION OF THE 
TIME-DEPENDENT STOKES EQUATION 

For a low Reynolds number fluid whose velocity and 
pressure are of the form of Eq. 0, the time-dependent 
Stokes equation is given by 



iwpv(q, z, uj)e lqx — — rj{uj)V 2 [v(g, z, ui)t 
+ V[P(q,z,w)e i « x ] , 



(Al) 



where p is the fluid density and rj(oj) is the frequency- 
dependent fluid viscosity. We also assume that the fluid 
is incompressible, 

V-v(a;,t)=0 => v x (q,z,u) = ~v' z (q,z,uj), (A2) 

where the prime indicates a derivative with respect to the 
z-coordinate. Using the identity V x [Vip(x, z, t)] = (for 
any scalar function tjj), we can take the curl of Eq. (|A1(I 
to eliminate the pressure, 

iujpn y (q,z,uj) +T)(w) [d 2 -q 2 ] n y (q,z,uj) = 0, (A3) 

where fl(x, z, t) = V x v(x, z, t) = fl y (x, z, t)y is the fluid 
vorticity, 

i 

n y (q,z,u>) = -v"(q,z,uj) - iqv z (q, z, u>). (A4) 

The solution to Eq. I|A3[1 can be written as 

£l y (q, z, uj) = 2qC\ cosh(/cz) + 2qC2 sinh(/cz), (A5) 
where C± and Ci are constants of integration and k = 



q 2 — ^jffi- The z-component of the velocity obeys the 
differential equation Eq. l|A4Jl . which has the solution 



v z (q, z, lo) —C^ cosh(qz) + C4 sinh(qz) (A-6) 

poo 

+ / dz'G(z - z')n y (q,z',uj), 



where the first two terms are the homogeneous solution 
to Eq. l|A4jl . and G(z) is the Green's function for the 
operator i(d 2 — q 2 )/q, G(z) = —ismh(qz)9(z), 9{z) being 
the Heaviside step function. Then the velocity is given 
by 

v z (q, z,lj) = 2C[q 2 cosh(kz) + 2C' 2 q 2 sinh(fcz) (A7) 
+ (C 3 - 2C[q 2 ) cosh(gz) + (C 4 - 2C' 2 kq) sinh(gz), 

where C] = Cj/(k 2 - q 2 ) for j = 1, 2. 



The pressure can be obtained from the x-component 
of Eq. PTJ) . 



P(q,z,u) 



[v' z (q,z,uj)-k 2 v' z (q,z,uj)]. (A8) 



Note that the z-component of Eq. (|A1|) is automatically 
satisfied by this pressure and the velocity field given by 
Eq. ETt . 

For overdamped fluids, where the inertial term in the 
time dependent Stokes equation - that is, the left-hand 
side of Eq. (|A1|I - is negligible, we can set p — 0, which 
causes k — > q. In this case, the velocity Eq. (|A7|) is given 
by 

V z (q, z,w) =C3Cosh(gz) + (C4 — C2)sinh(gz) (A9) 
+ C±qz sinh(qz) + C^gz cosh(qz). 



APPENDIX B: THE 
FLUCTUATION-DISSIPATION THEOREM 

The fluctuation-dissipation theorem relates correla- 
tions between the orthogonal dynamical variables de- 
scribing some system at equilibrium to the response of 
these variables to external forces. For the single fluid 
layer system we consider in this paper, there is only one 
dynamical variable, the height of the fluid, h(x,t). For 
the double layer system, there are two orthogonal dy- 
namical variables, the normal modes amplitudes h a (x, t), 
where a = ±. In both cases, the response function can be 
obtained from the normal stress boundary condition(s) at 
the fluid interface(s), which in the normal mode basis can 
be written in the form 



h a (q,uj) 



(Bl) 



The fluctuation-dissipation theorem relates the response 
function £ Q (q, w) to the height-height correlation function 



S a (q,t) = 



k B TL 2 



(h a (q,t)h* a (q,0)) 



(B2) 



(IlO 

2iriuj 



4 [Z a (q )( j)-&(-q,c 



where L 2 = j d 2 x and the star indicates complex conju- 
gation. 

The integral in Eq. can be evaluated by contour 
integration using a semi-circular, negatively oriented con- 
tour C whose arc lies in the negative imaginary half-plane 
and whose diameter is the real axis. We can see from 
Eq. I)B1|) that the frequencies for which ^~ 1 (q,uj) = 
are precisely the normal mode frequencies u> n [i.e. the 
frequencies that satisfy the normal stress boundary con- 
ditions Eq. (|Bip for f a = 0] . Because the system is over- 
damped, these frequencies lie on the negative imaginary 
axis, within the contour C; therefore, each normal mode 
frequency contributes a non-zero residue to the contour 
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integral. Furthermore, it can be shown that, for all sys- 
tems we consider in this paper, £*(— q,oS) — £ a (q, — w). 
This implies that the second term in Eq. I)B2|) has no 
poles inside the contour, nor does it contribute to the 
residues at the poles u) n . Assuming the poles at each co n 
are simple poles (which can be shown to be true in all 
cases), Eq. <|B2ll becomes 



S a (q,t) = lim 



[u - uj n ) e 



(B3) 



The time evolution of the surface height now depends on 
the entire deformation history of the surface convoluted 
with an exponential memory kernel. The simple form of 
the memory kernel is an artifact of the simplicity of the 
Maxwell fluid, but the general structure of this equation 
of motion holds for any more complex model of viscoelas- 
ticity. 

Taking a time derivative of Eq. ljC4fl we may write a 
single, second order differential equation for the ampli- 
tude of the height field at wave vector q: 



APPENDIX C: NUMBER OF NORMAL MODE 
DECAY RATES 

The number of decay rates in a dynamical system is 
simply related to the number of independent degrees of 
freedom |24|. It is therefore somewhat counterintuitive 
to observe two decay rates for capillary waves on a single 
supported viscoelastic layer and either four or five decay 
rates in the double layer system. It is perhaps most sur- 
prising to find the number of decay rates to change upon 
inverting the order of the layering of the two materials. 
In this appendix we explain this result in a simple way by 
considering the equations of motion in the time domain 
rather than in the frequency domain, as we have done 
throughout the rest of this article. 

Let us first consider the single layer case. Here and 
throughout this appendix we suppress any wave number 
q dependence for clarity. From the boundary condition 
Eq. I|1U|) we see that v z (z,uj) oc —iwh(u). Using the 
boundary conditions Eqs. (|5|)- (|llfl . we write the velocity 
in the form 



v z (z,u>) = —iu)h(uj)tp(z). 



Taken in combination with the stress boundary condition 
Eq. O we find 



^q 2 h(uj) = -iwr)(u)h(u)f(0), 



(C2) 



where f(z) = 3ip'(z) - ^"'{z)/q 2 . 

For a Newtonian fluid, r)(uS) — r\. Converting Eq. (|C2|) 
back to the time domain we recover a first order differ- 
ential equation for the decay of the surface height field: 



iq 2 h(t)=r)f(0)h(t). 



(C3) 



Thus, we expect only one decay rate at a given wave num- 
ber for capillary waves on a single supported Newtonian 
fluid. 

For the single supported Maxwell fluid the situation is 
more complicated. Now the viscosity is a complex fre- 
quency dependent function, which is given by Eq. (I23|l . 
As a result, conversion of Eq. 1|C2(1 into the time domain 
yields an integro-diffcrcntial equation: 



77fiU(t) = Voh(t) ~ — I dt'h(t')e 

7(0) T J_oo 



(C4) 



iq 
/(o) 



rh(t) + h(t) 



rh(t) + h(t) -E h(t). (C5) 



From this second order differential equation, which has 
two independent exponentially decaying solutions, we 
know to expect a double exponential decay for capillary 
waves on the supported Maxwell fluid. 

We now consider the case of two fluid layers. As in 
the single layer case we can use the boundary conditions 
Eqs. (|31fl - H3t)|) to write the fluid velocities in the form 

Vk,z(z, w) = -iw [hx(u))rp k i(z) + M^)^^)] , (C6) 

with k = 1,2 indexing the fluid layers. There remain 
three more boundary conditions. One enforces tangen- 
tial stress continuity at the fluid/fluid interface; the other 
two set the difference in the normal stresses at both the 
fluid/fluid interface and the free surface in terms of their 
respective surface tensions. These conditions may be 
written as, respectively, 

rn(u) J2 9ih(P)hk(u) = rn(w) J2 <?2fc(0)/ifcM, (C7) 



fc=l,2 



fc=l,2 



(CI) and 



11 q 2 h 1 (uj)= {-l) j iuh h {u) m {u)f jk {d 2 ) (C8) 

l2q 2 h 2 {oj) = ^ iLur) 2 (ui)g2k{0)hk(uj), (C9) 



fc=l,2 



where f jk (z) 



3Vk(s) " V*{*W and g jk (z) 



q 4>jk(z) + ip" k {z). In a manner analogous to the one 
used above with the single layer problem, we can rewrite 
Eqs. I|C7() - I|C9D in the time domain. Note that the pres- 
ence of the ioj factor in Eqs. (|CJ8I) and (|C9|I implies a time 
derivative. 

When both fluids are Newtonian Eq. I|C7|I involves no 
time derivatives, whereas both Eqs i|C8(l and i|C9(l are 
first order differential equations. Having a system of two 
first order differential equations, we expect two normal 
mode decay rates for a double layer of two immiscible 
Newtonian fluids. 

If either the upper or buried fluid is a Maxwell fluid, 
then it clear from the arguments given in the single layer 
case that Eq. I|C7|I becomes a first-order differential equa- 
tion in the time domain, whereas Eq. (|C8|I becomes a sec- 
ond order ODE. Eq. I|C9|) . on the other hand, depends 
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only on the viscosity of the upper fluid, 772, so it becomes 
a second order differential equation only when the up- 
per fluid is a Maxwell fluid; otherwise, it is a first order 
differential equation. Therefore, if the upper layer is a 
Newtonian fluid we have a set of two first order differen- 
tial equations and one second order differential equation; 



we expect there to be four decay rates in this case. If 
the upper fluid is a Maxwell fluid and the lower one is 
Newtonian, however, we now have a dynamical system 
described by one first order differential equation and two 
second order ones. In this case there will be five decay 
rates. 
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